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^ . Abstract 

o . 

■ Understanding how stochastic and non-hnear deterministic processes interact is a major challenge 
^Nj ' in population dynamics theory. After a short review, we introduce a stochastic individual-centered 

» I , particle model to describe the evolution in continuous time of a population with (continuous) age 

I' and trait structures. The individuals reproduce asexually, age, interact and die. The 'trait' is an 

' individual heritable property (d-dimensional vector) that may influence birth and death rates and 
interactions between individuals, and vary by mutation. In a large population limit, the random 

' process converges to the solution of a Gurtin-McCamy type PDE. We show that the random model 

' I has a long time behavior that differs from its deterministic limit. However, the results on the limiting 

^ PDE and large deviation techniques a la Freidlin-Wentzell provide estimates of the extinction time 

' and a better understanding of the long time behavior of the stochastic process. This has applications 

, to the theory of adaptive dynamics used in evolutionary biology. We present simulations for two 

• ■ biological problems involving life-history trait evolution when body size is plastic and individual 

' growth is taken into account. 
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a\ ■ 
o; 

00 ^ In so-called structured populations, individuals differ according to variables that influence their sur- 

' vival and reproduction abilities. These variables are morphological, physiological or behaviorial traits 

^ \ that can be inherited from parents to offspring and may vary due to genetic mutation. Other types 

■ of variables include ages, which are increasing functions of time, sex, or spatial location. This paper 
r> \ considers the dynamics of populations with trait and age structures. 

■ Age variables sum up individuals past histories and it may be biologically relevant to take several ages 
into account to model the present dynamics. Examples include the physical age (the time since birth), 
the biological age (the intrinsic physiological stage that can grow nonlinear ly in time), the age or stage 
of an illness, and the time since maturation. Taking age structure into account is required to study the 
evolution of individuals' life histories. For example, the salmon Oncorhynchus gorbuscha breeds when 
about 2 year-old, after returning to native freshwater streams, and most individuals die after this unique 
reproduction episode. This striking life history remains incompletely understood from an evolutionary 
point of viex, due in part to the lack of appropriate modelling tools. Moreover, age and trait structures 
interact as they change through time. This is because the age of expression of a trait and the life history 
shaped by the trait will influence the population dynamics, and hence the selection pressure on the trait, 
and its long term adaptive evolution. 

We start with a brief review of the literature on stochastic and deterministic models of populations 
with age structure in Section [TJ In Section [2 the links between stochastic and deterministic approaches 
are established: from a stochastic individual-based model (IBM), we recover the PDEs of demography 
via a large population limit. Section [2.41 proposes some illustrations of how deterministic and stochastic 
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approaches complement each other. The PDEs help study the measure-valued stochastic differential 
equations (SDEs) that arise from stochastic IBM. The SDEs in turn provide results that may not hold 
in large populations. An application to the study of extinction time is presented. Using these results, 
a second application to the so-called adaptive dynamics theory of biological evolution is developed in 
Section \TM This application highlights the importance of probabilistic modelling of point events. We 
address some interesting biological applications through simulations in Section [31 



1 Populations with continuous age (and trait) structure: a brief 
review 

In this brief review of the literature, we focus on population models with continuous age and time. For 
discrete-time models with age or stage classes, we refer the interested reader to [HI [201 El]- 

1.1 Deterministic models 

Population models with continuous age and time that generalize the equations of Malthus [62j and Ver- 
hulst's well-known logistic equation 87] have often been written as PDEs {e.g. [20 l l45 l [75 l [84 l [89]). 

PDE models for populations structured by a scalar age have been introduced by Sharpe and Lotka 
[5D] . Lotka [HI], McKendrick [H] and Von Foerster [37j: Vt > 0, Va > 0, 

— (a,t) + —(a,t)^—d(a)n(a,t), n(0,t)= / b(a)n(a,t)da, n(a,0) ~ naia), (1) 
ot da Jo 

where d, b and no are nonnegative continuous functions (death and birth rates, initial age distribution) 
with b bounded and no integrable. Equation H]) is known as McKendrick- Von Foerster PDE. It is a 
transport PDE that describes the aging phenomenon, with death terms and births on the boundary 
o = 0. A shortcoming of these models is that they do not take into account regulations of the population 
size by interference interactions or environmental limitations. 

The first nonlinear models for populations with age structure have been introduced by Gurtin and 
MacCamy |42]. The rates b{a) and d{a) of ([T]) become functions b{a,Nt) and d{a,Nt) of the scalar age 
a G M_|_ and population size Nt — Jq °° n{a, t)da G R+ at time t. Thus, the nonlinearity introduces a 
feed-back term that controls the population size. A particular case involves 

d(a, Nt) = d{a) + r^Nt (2) 

where d{a) is the natural death rate at age a and rjNt is the logistic competition term. The stationary 
solutions, their global stabilities and estimates of the rates of decay were studied by Marcati [63]. Stability 
conditions for a more general class of birth and death rates have been obtained by Farkas and 
Farkas [3S]. In [35], functions b(a,Nt) and d(a,Nt) are taken in C°'^(R^,R+) with the assumption: 
3d > 0, V(a, N) G IR+, d > d, and the characteristic equation whose roots determine the stability of the 
stationary solution is derived. In ^5], sufficient conditions for the stability of the stationary solutions are 
obtained for the cases where b{a,Nt) — bi{a)b2{Nt), in which the age-dependent fertility is weighted by 
the interactions with other individuals. 

The logistic competition term assumes that each competitor experiences the same interactions. Busen- 
berg and lannelli [TT] considered death rates that take the form d{a) + F{t,Nt, Si{t), . . . Sm{t)) where 
V« G lljTi] , Si{t) = J^°° ji{a)n{a,t)da, F is continuous in t and Lipschitz continuous with respect to 
the other variables, and the 7^ are positive continuous bounded functions that describe the competitive 
pressure exerted by an individual depending on its age. However the interaction term remains non-local 
(and thus finite-dimensional). Local interactions mean that the competition terms depend on integrals of 
the form J^°° U {a, a)n{a, t)da where U{a,a) is the interaction exerted by an individual of age a on an 
individual of age a. This case is covered by the work of Webb [52], who proposes PDEs with coefficients of 
the form b{a, n(., t)) and d{a, n{., t)), where n{., t) belongs to the Banach space L^(R+, R+) of integrable 
functions w.r.t. a G R+. 

More recently, entropy methods have extended the notion of relative entropy to equations that are 
not conservation laws. This approach (Michel et al. [7T], Mischler et al. [72], Perthame [75 , Perthame 
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and Ryzhik [76] ) provided new proofs of existence of solutions, and new insights into long time behavior 
for instance. 

There are many examples in the literature where additional structures other than age-based have 
been taken into account. The population can be divided into a finite number of classes: see Webb [89] for 
a survey on multi-type models, including prey-predator and epidemiological models (e.g. [T^ [55] ). 
In epidemics models, considering the classes of susceptible, infectious and/or removed individuals with 
different interactions between these groups, leads to different propagation dynamics of the infection. Age- 
structure is important as infection or detection rate may depend for instance on the time since the first 
infection. 

To our knowledge, populations with age structure and spatial diffusions (where motion is described 
with a second order differential operator) were first studied by Langlais [5S1 EH US] and Busenberg 
and lannelli [10] . Structure variables with values in a non- finite state-space other than position have 
been considered by Rotenberg [75] for instance (for multitype populations, see e.g. [SI]). In [75], the 
maturation velocity a; G A" =]0, 1] is defined as a trait that can change at birth and during the life of an 
individual. For an individual born at time c S M and with traits Xi €]0, 1] on [ti, {i E N, to = c and 
(ii)igN are times of trait change in [c, -l-oo[), the biological age at t > c is a = X^ieN ^ii'^i+i A t — A t). 
Here, the biological age a differs from the physical age t — c as soon as 3i G N, ti < t and Xi 1. A 
nonlinear version of this model has been studied by Mischler et al. ^72, . 




1.2 Stochastic models 

Age structured branching processes that generalize the Galton- Watson process [41] have been studied by 
Bellman and Harris [4l|43], and then Athreya and Ney ([2J Chapter IV). In these non Markovian models, 
the lifespan of an individual does not follow an exponential law. Upon death, a particle is replaced by a 
random number of daughter particles, with a reproduction law that does not depend on the age of the 
mother, nor on the state of the population. 

The assumptions of births at the parent's death and of independence between the reproduction law 
and the age of the parent are biologically restrictive. Kendall f54^. Crump and Mode [221 [23], Jagers 
[151 H5] , Doney [31] , study birth and death processes in which an individual can give birth at several 
random times during its life, with rates that may depend on its age. These processes are called age- 
structured birth and death processes. More recently Lambert [57] used contour processes to study the 
properties of splitting trees which are formed by individuals with independent and identically distributed 
(i.i.d.) lifespans and who give birth at the same constant rate during their lives. Lambert's results 
provide a new interpretation of the link between branching processes and Levy processes |40[ 177] . 

In the two types of models introduced above, individuals alive at the same time are independent, which 
is also a biologically restrictive assumption. Wang [SH], Solomon [81 have considered birth and death pro- 
cesses in which the lifespans of individuals are independent, but in which the birth rate and reproduction 
law of each individual depend on the state of the population. Oelschlager [73 , Jagers [48 , Jagers and Kle- 
baner S2j generalize these models by including interactions in both the birth and death rates, but these 
rates remain bounded. In the models mentioned in this paragraph, the population at time t is discrete 
and represented by a point measure of the form Zt{da) = X^iL'i ^ai{t){dct). Each individual is described 
by a Dirac mass charging its age ai(t), and Nt is the size of the population. The birth and death rates con- 
sidered by |3H1I521[7H] have the form b{a, {Zt, U{a, .))) and d{a, {Zt,V{a, .))), where U and V are bounded 
interaction kernels, b and d are bounded functions, and we use the notation {Z, f) — /q^°° f{a)Zt{da). In 

[52] 148] for instance, the interactions vanish in the limit: lim^^ b(a, /q^°° U{a, a)Z{da)) = 6(a) and 

the same is true for the death term, so that the authors recover population behaviors that are independent 
of population size in the limit of large populations. These models thus exclude logistic interactions. Such 
interactions are included in the models considered in [M] . 

Interactions between individuals imply that the key branching property ceased to hold, and classical 
approaches based on generating or Laplace functions {e.g. |2]) do not apply anymore. However, in a 
large population limit, a law of large numbers can be established |HH1 IHIl [731 [HS] whereby a limiting 
deterministic process identified as a weak solution of a generalization of the McKendrick-Von Foerster 
PDE ll]) is obtained. We present a synthesis of these results in Section [2] Once the link between PDEs 
and stochastic processes is established, results for PDEs prove to have interesting implications for the 
stochastic process, even if the behaviors of the stochastic model and of its corresponding deterministic 
limit differ (see Sections 12.41 and l2.5|) . 
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Let us mention that there are other large-population scalings that have been the focus of abundant 
research. Superprocess limits of age-structured branching processes (where the mother gives birth at 
her death) are obtained by rescaling the lifespan and mass of the particles, by increasing their number 
and by modifications of the birth rate or reproduction law {e.g. Dynkin [33], Kaj et Sagitov [53|). For 
age-structured birth and death processes (where the mother gives birth at random times during her life) , 
similar approaches were taken by Dawson et al. [25 and Bose and Kaj • 

Multitype populations with age structure have been considered by many authors {e.g. Athreya and 
Ney [3], Chap. V. 10). A stochastic model of epidemics with age structure has been considered for instance 
in Clemengon et al. [21]. Stochastic models of populations with age and trait structures have been 
addressed by Jagers [SUl |5T] and Meleard and Tran [55] for instance. In [55], the mutation rate that 
is responsible for generating variation in the trait is decreased so that the mutations become more and 
more rare. As a consequence, the timescales of mutations and demography (births and deaths) become 
separate, which is a typical assumption of the evolutionary theory of adaptive dynamics (see Section f2.5p . 

2 From stochastic individual-based models of age-structured pop- 
ulations to deterministic PDEs 

2.1 A microscopic model 

We study a microscopic stochastic model, for which the dynamics is specified at the level of individuals. 
The model takes into account trait and age dependence of birth and death rates, together with interactions 
between individuals. This process generalizes the approach of Fournier and Meleard [38j . Champagnat 
et al. [THIEH]. When the population is large, we establish a macroscopic approximation of the process, 
which describes the evolution at the scale of the population (individual paths are lost). These results are 
taken from |86j and allow to link the probabilistic point of view to the deterministic approach reviewed 
in Section 1. 

The IBM represents the discrete population at time t > by a point measure 

= (3) 

1=1 

where each individual is described by a Dirac mass (x,a) e A" := A" X K+, A" C R"* being the trait 
space and M+ the age space. We denote by Nt the number of individuals alive at time t. The space of 
point measures on X is denoted by A4p{X) and embedded with the topology of weak convergence. For 
a function /, we denote by {Zt, /) the integral f{x, a)Zt{dx, da) = X^iJi f{^ii'^)i An individual 

of trait X and age a in a population Z G A4p{X) reproduces asexually, ages and dies. When reproduction 
occurs, the trait is transmitted to offspring unless a mutation occurs. The mechanism is modelled as 
follows: 

• Birth rate is b{x, a) e R+. With probability p e [0, 1], the new individual is a mutant with trait x+h 
where h is chosen in the probability distribution k{x,a,h)dh (dependent upon the characteristics 
of the parent) where dh is the Lebesgue measure, 

• Death rate is d{x,a, ZU{x,a)) E M+. It depends on the trait and age of the individual but also 
on the measure Z describing the whole population. The function U : i— s- R™ is an interaction 
kernel: each of the m component of C/((a;, a), (y, a)) describes a different interaction of {y^a) on 
(a;, a), and ZU{x, a) — U{{x, a), (y, a))Z{dy, da). 

• Here, aging velocity is 1. 

Assumption (HI) The birth rate b{x, a) is assumed to be continuous and bounded by a positive constant 
b. The function k{x, a, h) is assumed to be bounded by a positive constant k. The death rate is assumed 
to be continuous in (a:, a), Lipschitz continuous in the interaction term, lower bounded by a strictly 
positive constant and upper bounded by d{l -\- {Z, 1)) with d > 0. 

Remark 1. 1. More general models are considered in [551 [HS], where vectorial ages with nonlinear aging 
velocities are allowed. 
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2. For appropriate choices of b and d, some of the raodels described in Section FOl can be recovered. If 
b{x, a) = 6(a) and d{x, a) = d{a), we recover the age-dependent birth and death processes of [12 [231 EH 
[48l[54]. The choice b{x,a) — b[a) and d{x, a, ZU{x, a)) = d{a) + rj{Z,l) gives the rates used by [42] , 
Considering bounded functions d and vectorial functions U leads to models similar to those considered 
in [51113]. 

Following [3H1 [13 [m [HS] , we describe the evolution of (.^t)tgB+ by a SDE driven by a Poisson point 
process. The mathematical difficulty lies in the fact that the birth and death rates depend on time 
through age. ^ 

Let Zq g Mp[X) be a random variable (r.v.) such that E((Zo,l)) < +00, and let Q{ds,di,d9,dh) 
be a Poisson point measure (P.P.M.) on M+ x £ := M+ x N* x x X with intensity q{ds, di, dO, dh) := 
ds®n{di)®d9® dh and independent of Zq (for the general theory of P.P.Ms, see for instance [17] Chapter 
I). 

Let us denote by Xi{t) and Ai{t) the trait and age of the individual at time t, the individuals 
being ranked in the lexicographical order on R'' X R , (see |38[ 186] for a rigorous notation) . 



i=i 



^(Xi(s_),t-s)l{0<e<mi(s,Zs_ ,i,h)} 



+^(Xi(s_)+/i,t-s)l{mi(s,Zs_ A,h)<e<m2(s,Z,_ ,i,h)} 

— ^{Xi(s-),Ai(s-)+t-s)'^{m2{s,Z^_,i,h)<e<m3{s,Z,_ ,i,h)} 



Q(ds,di,d9,dh), (4) 



where: 



mi{s,Zs_,i,h) = {l-p)b{X,{s^),A{s^))k{X,{s^),A{s^),h) 

W2(s, Zs_ ,i,h) = TOi(s, Zs_ , i, h) + pb{Xi{s^), Ai{s^))k{Xi{s-), Ai{s-),h) 

mais, ,i,h) = 1712(3, Z,_ ,i,h) + d{X,{s-), A,{s-), Z,_U{X^{s-), A,{s-)))k{X^{s-), A,{s-), h). 

The interpretation is as follows. To describe the population at any given time i, we start with the 
individuals present at t = 0. Then we add all the ones that were born between and t and finally we 
delete the Dirac masses corresponding to individuals who died between and t. The individuals of the 
initial condition have at time t an age increased by t (first term of ([3])). When a birth takes place at time 
s, in the measure describing the population at time t we add a Dirac mass at age t — s (the age expected 
for this newborn individual, second and third term of ([1])). When an individual of age a dies at time s, 
in the measure describing the population at time t we suppress a Dirac mass at age a -\- 1 — s (the age 
this individual would have had if it had survived). As the rates vary with time, we use a P.P.M. with 
an intensity that upper-bounds the age-dependent rates and use an acceptance-rejection procedure (the 
indicators in 9) to recover the rates that we need. 

We emphasize that exact simulation of this SDE can be done easily (see Section [3]) . The algorithm 
formally accounts for the computer programs that many biologists use to run simulations of individual- 
based population models {e.g. [25]). 

Proposition 1. Under (HI ), for every given Poisson point measure Q on R_|- x £ with intensity measure 
q, and every initial condition Zq such that E((Zo,l)) < -t-00, there exists a unique strong solution to 
SDE in D(R+, A^p(A:')). The solution is a Markov process with infinitesimal generator given for all 
Ff{Z) ^ F{{Z,f)) with f e C°'^{X,R), F e Ci(M,R), and Z e Mp{X) by 

LFfiZ) = / [da fix, a)F' {{Z, /)) + {Ff {Z + 5(,,o)) - {Z)) b(x, a)(l - p) 
Jx 

+ / {^f + kx+h.o)) - Ff {Z)) b{x, a)p k{x, a, h)dh 

+ {Ff {Z - 5(,,,)) - Ff {Z)) d{x, a, ZU{x, a))] Z{dx, da). (5) 

For the proof, see Propositions 2.2.5, 2.2.6 and Theorem 2.2.8 in 85j. The result is based on moment 
estimates and uses the algorithmic construction of ([4]) that will be detailed in Section [3] The behaviour 
of the solution is far more difficult to study. For this reason, and also in order to link the SDE ^ with 
the PDEs that are classically introduced in demography (see Section [TTT]) . we consider a large population 
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limit. Thus we let the size of the initial population grow to infinity proportionally to an integer parameter 
n while individuals and interaction intensities are assigned a weight 1/n. Under proper assumptions on 
the initial condition, the sequence of rescaled processes (Z"')„gN converges, when n +00, to the solution 
of a deterministic equation identified as the weak form of a PDE ; this PDE generalizes the deterministic 
models described in Section [TTTl 



2.2 Large population limit 

The large population renormalization that we consider is inspired by the work of Fournier and Meleard 
[55] . Full details and proofs can be found in [Ml \EE\- We let the size of the initial population tend 
to infinity proportionally to the parameter rt £ N* and renormalize the weights of the individuals and 
their interaction by 1/n. Following on Metz et al. [70 , n was named "system size" by Champagnat 
et al. |19J. More precisely, we consider a sequence (Z", t e R+)„gN* such that nZ^ satisfies SDK 
(HI) with J7" = U/n so that the interaction term ZU{x,a) = ^((^' (^i' ^j)) replaced with 

i X^i^i U{{x^ a), [xi^ a-i)) ~ Z'^U{x^ a) and such that: 

Assumption (H2): The sequence (Zq)^^* = ^(^i(o),Ai(o))) converges in probabihty to the 

measure ^0 belonging to the space A^f('^) of finite measures embedded with the weak convergence topol- 
ogy- 

This large population rescaling can be understood as reflecting resource limitation as the size of the 
population increases: we have to assume that the biomass of individuals decreases if the system is to 
remain viable. As the individuals become smaller, their interactions decrease proportionally to their 
mass. 

Proposition 2. Under (HI) and (H2), 

(i) //sup„gN. E((Z5', 1)2) < +00, then for every f e Ci'i'°(M+ x X ,R) the process 

M,"^/=(Zr,/(.,.,i))-(Zo",/(.,.,0))- f f [daf{x,a,s)+dJ{x,a,s) + f{x,0,s)bix,a){l-p) 

Jo JX 

+ / f{x + h,0,s)b{x,a)pk{x,a,h)dh — f{x,a,s)d{x,a,Z^U{x,a)) Zl'' {dx , da) ds , (6) 

is a square integrable martingale with the following quadratic variation process: 

/2(x, 0, s)6(x, o)(l — p) + / f^{x + h,0, s)b{x,a)pk{x,a,h)dh 
+ f^{x, a, s)d(x, a, Z^U{x, a))] Zl'{dx, da) ds. (7) 

(a) If 3ri > 0, sup„gj^, E((Zq , 1)'^+'') < +00, then the sequence (Z")„gN* converges in distribution in 
D(M+, A^i?(A')) to the unique solution ^ G C(M+, A^i?(A')) of the following equation: V/ : {s,x,a) i— > 
fsix,a) e Ci^O'i(R+ X ^^^), 

(6, /(•, i)> =(^0, /(., .,0)) + / / [dafix, a, s) + d,f{x, a, s) + f{x, 0, s)b{x, a){l ~ p) 



1 

It = - 
n 



Jx 



10 JX 

+ I f{x + h,0,s)b{x,a)pk{x,a,h)dh — f{x,a,s)d{x,a,£^sU{x,a)) 



^sidx,da)ds, (8) 

This result is proved by a tightness- uniqueness argument (see [3H1 \SE\)- Point (i) is obtained by 
stochastic calculus for jump processes. Heuristically, since the quadratic variation is of the order of 1/n, 
it vanishes as n ^ -l-oo. Moreover, when births or deaths occur, we add or delete individuals of weight 
1/n. This explains why in the limit we obtain a continuous deterministic process. 

The link between ^ and PDEs is specified in the next proposition (for proofs see Prop. 3. 6 and 3.7 
in [HS]). 

Assumption (H3): ^0 admits a density no(a;, a) with respect to dx da on X. 
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Proposition 3. Under (HI), (H2) and (H3), the measures S M.f{X), fort G M+, admit densities 
n{x,a,t) with respect to dx®da on X. The family of these densities .,t))teR+ a weak solution of: 



dtn{x,a,t) = —dan{x,a,t)~dyx,a,J__U{{x,a),{y,a))n{y,a,t)dydajn{x,a,t) (9) 

n{x,0,t) = / n{x,a,t)b{x,a){l — p)da + / b{x — h,a)pk{x — h,a,h)n{x — h,a,t)dhda 

n{x,a,0) = no(x, a). 

These equations generalize tlie McKendrick-Von Foerster's PDEs. Tliey describe ecological dynamics 
at the scale of the population (individual trajectories are lost). The densities n{x, a, t), when they exist, 
correspond to the number density in the sense of Desvillettes et al. |27j , Champagnat et al. jl8| , which 
describes the trait and age distributions of a "continuum" of individuals. Here, existence and uniqueness 
of a weak solution of ^ are obtained by probabilistic proofs. An interesting property of equation ([5]) is 
that function solutions do not always exist (see [86J). 



2.3 Central limit theorem 

In large populations, the microscopic process can be approximated by the solution of PDE ([9]) that 
generalizes classical PDEs of demography for populations in continuous time and structured by a scalar 
age. In order to construct confidence intervals or to assess the quality of the approximation, it is useful 
to study the fluctuation process, 

Vte [0,T],VneN*, T]];-{dx,da) ^y^{Zl'{dx,da) - ^t{dx,da)) . (10) 

For n g N*, 77" is a process of ©([0, T], A4s{'^))- Since the space Aisi^) of signed measures on X cannot 
be meterized when embedded with the topology of weak convergence, we follow the works of Mctivier 
[55] and Meleard [H^, and consider 77" as a distribution- valued process. 

Let us denote by the space of smooth functions with compact support. For a multi-index k e N'^^^, 
we denote by f the derivative d^'^^f/dx^'^ ■ ■ ■ dx^/ da^'^+^ . The Banach space C^'''' is the space of 
functions / of class such that for k such that \k\ < /3, lim D''f{x, a) / (l+a^) = when a +00 or when 
X — > dX, embedded with the norm H/Hc^.t = ^\k\<i3^^P{x a)ex +^^)- The Hilbert space 

Wq''' is the closure of with respect to the norm |1/|1^;3,^ = Ix^\k\<i3 \D''f{x, a)p/(l + a^^) dxda, 

being the space of smooth functions with compact support. We denote by C~^''^ and Wq^'^ their 
dual spaces. We are going to consider the following embeddings. Let Z) = [(d + l)/2] + 1, 

^3D+1.0 ^ VFo'^+l,^ ^ VK2^+1^2I3 ^ ^D+l,2D ^ (jD,2D ^ ^D,3D ^ ^0.313+1 

the second embedding being Hilbert-Schmidt (e.g. 1\ p. 173). 

The tightness is proved in the space Wq^^^^'^ using a tightness criterion for Hilbert-valued processes 
that is stated in [H7| (Lemma C), and which requires the Hilbert-Schmidt embedding Wq^'^^'^ ^h.s. 
Wq^~^^'^^ . (7"^^+^-''^^ and C~^'^^ are the spaces in which the estimates of the norm of the finite varia- 
tion part of the process are established (see Section 4.4.2 of 85J). We need these two spaces since da maps 

norm of the martingale part is controlled in Wq Finally, uniqueness of 

the limiting value is proved thanks to the embedding W(^^^~^^'^ ^ (j-3D+i,o^ 

Assumption (H4): We assume that Assumptions (HI), (H2) and (H3) are satisfied and that: 

• The sequence (??o )neN* converges in W^^'^^ to 770 and sup„£pj* IE(||f7o II^-d.3d) < +00. 

• sup„gN.E((/_^|a|6-Dzj(da;,rfa))^) < +00. 

• The functions (x, a) h(x, a), (x, a) i— > k{x, a, x'), {x, a) U {{x' , a), (x, a)) and {x, a) ^ U {{x, a), (x', a)) 
belong to for almost every {x',a) E X. 

• The death rate (x, a, u) 1-^ d{x, a, u) is such that for all u, it belongs to (72^+2,0 with a norm bounded 
by a polynomial in u. We assume that its derivative with respect to u is Lipschitz continuous. 
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Proposition 4. Let T > 0. Under Assumptions (H4), the sequence (7?"')neN* converges in distribution 
in B{[0,T],Wq^°'^^'°) to the unique continuous solution of: Vi S [0,r], V/ S 1^0^+^'-°, 



imJ) ^{voJ) + Wt{f) + 

Jo Jx 

d{x,a,^sU{x,a))f{x,a) 



da fix, a) + b{x, a) ((1 - p)f{x, 0)+p f{x + h, 0)k{x, a, h)dh) 



f{y, a)duiy, a, ^sU{y, a))U {{y, a), {x, a))^s {dy, da) 



X 



rjs{dx, da) ds 
(12) 

where (Wt{f))teWL+ « continuous centered square-integrahle Gaussian process with quadratic variation: 



{W{f))t = 



JX 



b{x,a){{l~p)f^{x,a)+p / f^{x + h,0)k{x,a,h)dh) 



X 



f'^{x, a)d{x, a, £^sU{x, a))] £,s{dx, da) ds. 



(13) 



Notice that a benefit from IBMs is that the stochastic approach with point measures suits the for- 
malism of statistical methods which involve sets of individual data {e.g. [3Tl[S]). This makes it possible 
to use the battery of statistical methods which deal with problems that may not always be treated with 
deterministic methods (such as missing or noisy data). The convergence and fluctuations of Proposi- 
tions [5] and [5] would then provide consistence and asymptotic normality for the estimators. Calibrating 
the parameters of a PDE thanks to its underlying microscopic interpretation provides an alternative to 
deterministic approaches {e.g. [SIE])- 



2.4 Extinction in logistic age structured populations with constant trait 

Extinction is one of the recurrent and important issue when studying the ecology of a population. Here, 
the conclusions given by the deterministic and the stochastic models differ, but we will see that they still 
provide complementary information. 

For the sake of simplicity, we consider here a logistic age-structured population without trait variation. 
We will write for instance no (a) instead of no{x,a) of Assumption (H2). An individual of age a in a 
population of size N gives birth with rate b{a) and dies with the rate d{a)+riN as in From Proposition 
[21 under the large population renormalization of Section \2.2\ the sequence {Z")neti* converges to the 
unique weak solution of the Gurtin McCamy PDE with the death rate ([2]). Let us recall some known 
facts about this PDE. We denote by n(ai, 02) ~ exp ( J^^ d{a)daj the probability of survival from age 
ai to age 02 in absence of competition. 

Proposition 5. Under (HI) and (H3), there exists a unique classical solution of class to the Gurtin- 
McCamy equation with death rate 0). It is given by: Va G IR+, Vi G ]R-|_ 

f ,x Nov{a,t) J no{a-t)n{a-t,a)/No if a>t 

n(a,t) = ; — x:— with via, t) — < \„, J .j. ^ (14) 

l + NoJl^J^^via,s)dads ^ ' \ B{t-a)U{a,0) if a<t ^ > 

where uq is the density of the initial condition, Nq its L^-norm and: 

B(t) = So * ( V.g*")(t) with g{a) ^b{a)Ii{Q,a)ta>a and Bo{t) = lt>o / b{a + t)^^^^^U{a, a + t)da. 
n=a ■Jo ^0 

The quantity B{t) which appears in the computations can be interpreted heuristically as the number 
of births at time t. For a proof, see Webb [HH] (Section 5.4). 

Since the family {n{a, t)da)t^M^ defines a weak solution of the Gurtin-McCamy PDE in this case, we 
obtain the following corollary of the uniqueness property stated in Proposition [2] and of Proposition [5l 

Corollary 1. In the case of logistic age- structured population, the limiting process ^ of the sequence 
(Z")„gN* satisfies: G M+, S,t{da) — n{a,t)da where n{a,t) is explicitly given in 



We have used here the results known for the PDE to obtain an explicit expression and regularities of 
the density of ^ obtained by probabilistic proofs (Propositions [2] and [S]) . 

As explained in Section [H the long time behavior of the Gurtin-McCamy PDE is well known. In 
particular: 



Proposition 6. (see \6'A \8 9l ) Under the assumption Rq :— J^°° b{a)ll{0, a)da > 1, which expresses 
the renewal of generations in absence of competition, there exists a unique nontrivial asymptotically 
exponentially stable steady state given by: 



fi{a) = -J— where Ai satisfies 1 = / e ^°6(a)n(0, a)c?a. (15) 

rj /q e-^i"n(0, a)da Jo 



100 200 300 400 500 1 00 200 300 400 500 1 00 200 300 400 500 

Time Time Time 

(a) (b) (c) 

Figure 1: Simulated size of a logistic age- structured population. The model is given by ^ rescaled with n = 1, 5 
and 10 (from left to right). 

Even if simulations (shown in Figure [1]) may lead us to the false impression that the long time be- 
havior of the stochastic processes is the same as for their deterministic limit, the situation is much more 
complicated. We show that the stochastic process goes extinct almost surely, but after having spent 
an exponentially long time in the neighborhood of the stationary solution ^{da) = n{a)da (|15|1 of its 
deterministic approximation. Here, deterministic and stochastic models yield different conclusions, but 
the results for the limiting PDF make it possible to understand the behavior of the stochastic processes. 
With small system size the population size demonstrates wild fluctuations, and extinction is likely to 
occur quickly. As system size increases, fluctuations in population size are reduced, extinction takes 
more time, and the mean size of non-extinct populations converges toward the equilibrium solution of 
the deterministic model approximation. 




First of all, on compact time intervals, the stochastic process belongs with large probability to a tube 
centered on the deterministic limit. This implies by Proposition [2] that: 

Ve > 0, 3te > 0, > ie, Hm P {p{Z'^, > = 0, (16) 

where p is for instance the Dudley metric which meterizes the topology of weak convergence on M.p{X) 
{e.g. [75] p79). Heuristically, stochastic and deterministic processes have the same behavior on compact 
time intervals. In particular, since the deterministic process enters a given e-neighborhood of its nontrivial 
stationary solution ^ after a sufficiently large time that does not depend on n, the stochastic process enters 
an 2e- neighborhood of ^ with a probability that tends to 1 when n — s- +oo. However, when T — s- +oo, 
the long time behavior differ. 

Proposition 7. For fixed n, one has almost sure extinction: 

P(3t e M+, (Z,",l) =0) = 1. (17) 

Because of logistic competition, the population stays finite with a size that is controlled and cannot 
grow to infinity (see Proposition 5.6 [86 ). The stochastic process finally leaves the neighborhood of ^ 
to drive the population to extinction. We can thus ask how long the stochastic process stays in the 
neighborhood of ^. 

Let < 7 < 1) and assume that we start our microscopic process in the neighborhood of the 
equilibrium ^: p{Zq,^) < 7 a.s. An estimate of the time T" = inf{t e R+, p{Z",£^) > 7} that the 
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stochastic process spends in the neighborhood of ^ can be estabhshed by large deviation results a la 
Freidlin and Ventzell [Ml [Mj and generahzed to our measure- valued setting (see Section 5.2 [HS]). The 
proof (see [86]) is based on the fine properties of the convergence of Z" to its deterministic limit. 

Proposition 8. 3V >Q,V_> 0, V(5 > 0, 

lim P(e"(^-*') < T" < e"(^+'')) = 1 (18) 

n — > + oo 

When n — > +00, T" tends to infinity in probability and extinction eventually disappears in the 
limit: unlike the microscopic process that gets extinct almost surely, its deterministic limit admits a 
nontrivial equilibrium. As a corollary, for every sequence (tn)„gN* such that lim„_>+oo 
and lim„^+oo tn — +00, 

hm V{p{Zl]J)>e)^Q. (19) 

Equation p8p is interesting in itself to understand the persistence time of a stochastic finite population 
; it is useful also to separate timescales in models of adaptive dynamics, as we shall see in the next 
sub-section. 



2.5 Application to the adaptive dynamics theory of phenotype evolution 

In the theory of phenotype evolution, mutations are responsible for generating trait diversity in the 
population while natural selection acting through interaction and competition determines the traits that 
become fixed. Important questions in evolutionary biology require to include age structure in population 
models, as Charlesworth ^20j, Medawar ^66^ and Stearns [H2j emphasized. The selective pressure and 
fixation probability of a trait may be functions of the age at which the trait is expressed, as well as of 
the age structure of the population. Reciprocally, the traits that are fixed can modify the age structure 
of the population. 

The so-called "adaptive dynamics" framework considers large population limits and rare mutations. 
Adaptive dynamics theory describes the evolution of traits in the population when it is possible to 
separate the timescales of ecology (birth and death events) and evolution (generation of new traits). 
For populations with only trait structure, adaptive dynamics theory has been introduced and developed 
by Hofbauer and Sigmund [33] , Marrow et al. [53] , Dieckmann and Law [5S] , Metz et al. [7D1 [32] and 
Champagnat et al. [TTllTOlfTH] . 

These studies have been generalized to populations with age and trait structures by Meleard and 
Tran ^68j . Therein, equations of adaptive dynamics are derived from microscopic models by following 
Champagnat [17j and several examples are studied and illustrated with numerical simulations. When 
mutations are sufficiently rare, natural selection wipes out the weakest competitors so that the resident 
population at the next mutation is monomorphic and in a stationary state. The estimates of extinction 
time provided in Section f2.4l tell us that mutations should not be too rare if we wish to neglect the events 
of extinctions before the occurrence of a new mutant. By neglecting the transition periods, it is possible 
to describe evolution as a jump process by considering only the successive monomorphic equilibria (trait 
value and corresponding stationary age structure). If there exists a unique nontrivial stable stationary 
solution S^x{da) — n{x,a)da, it is sufficient to consider the trait valued process that jumps from one 
equilibrium trait to another. This process, called trait substitution sequence process for age- structured 
populations (TSSPASP), introduced and analyzed in [68], can be described as follows. 

• In a monomorphic population of trait x at equilibrium, mutants are generated with rate 

pN^ 



p I b{x, a)n{x, a)da = pn{x,0) 



E(T,) 



where = /q n(x,a)da is the size of the population at equilibrium and E{Tx) is the expected 
individual lifespan: 



The probability that the resulting mutant population with trait x + h replaces the resident popu- 
lation with trait x cannot always be calculated explicitly when we deal with populations that are 
structured by age. However, it can be computed numerically {e.g. [68]). 
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The theory of random point processes is necessary to model these events that occur at discrete random 
times. Evolutionary models have also been developed from a deterministic point of view, but with different 
renormalizations that do not describe the same rare mutations (Diekmann et al. [30] , Carrillo et al. [14] ) . 

If mutation steps are additionally assumed small, the TSSPASP can be approximated by the solution 
of an ordinary differential equation that generalizes the " canonical equation of adaptive dynamics" (see 
Dieckmann and Law [33], Champagnat [TO]). 

Recently, models of populations with age structure have been considered in adaptive dynamics theory 
deahng with functional traits {e.g. [551 (Ml El])- The traits a;(a) are functions of age (growth curve, 
flowering intensity, birth or death curves...) In Dieckmann et al. [28] or Parvinen et al. [74 , the birth 
and death rates are functions of the functional trait, but not of age (for example, the averaged form 
b{x) — B{J^°° x{a)da) with B G Ch(M,R+) was assumed). 

3 Age and trait structures: model examples and insights from 
simulations 

We now present several examples of interesting biological questions involving populations with age and 
trait structures. Mechanisms are described at the individual level, and we are interested in the resulting 
macroscopic dynamics. Not surprisingly, the more realistic and complex the model is, the more diffi- 
cult it is to carry out an analytical study of the model behavior. Simulations offer an alternative and 
complementary approach. 



3.1 Simulation algorithm 

Following the algorithms proposed by Fournier and Meleard [35], it is possible to simulate exactly the 
law of the process Z ([3]), without approximation scheme or grid. The main difficulty comes from the fact 
that the rates depend on age that changes continuously in time. To handle this, we use acceptance-reject 
procedures. 

Let Zq E A4p{X) be an initial condition. We simulate by recursion a succession of birth and death 
events that modify the size of the population. Let us set Tq = for the event number 0. Assume that 
we have already simulated k events (fc G N), and that the last of these events had occurred at time Tk. 
The size of the population Nt^ = {Zt^. , 1) at time Tk is finite and the global jump rate is upper bounded 
by 6iVTfc + d(l + Nt^)Nt^, which is finite. To obtain Tk+i for the k + 1^^ event, we simulate candidate 
events from time T^. The latter are given by a sequence {Tk/)££ti of possible event times following a 
Poisson point process with intensity BNt^ -f- c?(1 -f NT^)NTf.. The first of these times which is accepted 
by the procedure defines Tk+i. 

0. We set Tk,Q Tk, Nr^ := Nt, and £ 0. 

1. We simulate independent exponential variables Sk e with parameter 1 and we define t^^ — Tk i + 

ek,,/[bNT,+d{l + NT,)NT,]. 

2. On the interval [rk.i, Tkj^+i [, only aging takes place. For i g |1, iVT^], the age of individual i becomes 
Ai{Tk,i) + Tkj+i - Tkj at time Tk^i+i- 

3. We simulate an integer valued r.v. Ikj uniformly distributed on |1, Nt^\, and we define the following 
quantities in [0, 1]: 

~ , r \ (1 -p)K^^fc.^(-f(fc,« + l)-)>^Ife,£(T-(fc,« + l)-)) 

^^iW,Zr,,,,^,,_,h,i) - h + dil + Nr,) 

m2(Tk,e+i,Zr^^_^^-^^^,lk,e) = mi{Tkj+i, Zr^i^i.^-^^_ , hj) H b + ij(l + jVy,) 

d (-''^-rfc,£('r(fc,<:+i)-)>^i"fc,(!('^(fc,<!+i)_)i-^T(fc,^+i)_ f/(-'f/fc,£(-r(fc,^+i)_), A7^_^(r(fe,f+i)_))^ 



h + d{l + Nt, 



4. We simulate a r.v. Qk.e with miiform law on [0, 1 
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(a) If < Qk,i < 'miiTk,e+i, Zt^^. , Ik,e), then the k + 1 event occurs: Individual Ik/ gives 
birth to a clone with age G and one defines Tk+i — Tkj+i- 

(b) l^■ml{Tk,l>+l,Zr^^ ^^^^_,Ik,l) < Qk,i < m2(rfc,£+i, , 4,^), then the /c + 1*^ event occurs: 
Individual Ik,e gives birth to a mutant of age and trait x + h where h follows the distribution 
probability of density k{Xi^_^{T(^k,t+i)~), Ai^ ^{T(^k/+i)-), h'). One defines Tk+i = Tkj+i- 

(c) lim2{Tki+i,Zr^^ ,^^^_,Ik,i) < Qki < ma (rfe,^+i, , 4^^), then the /c + l*^ event occurs: 
Individual Ik/ dies and one defines Tk+i — Tk.i+i- 

(d) If m^{Tk.t,+i, ^T(fc (.+!)_ J Ik.i) < ©fe/ then nothing happens. We reiterate the algorithm from 1. 

th 

with ^ + 1 in place of £ until we obtain the k + 1 event. 

Remark 2. To each individual are associated three clocks corresponding respectively to birth (with or 
without mutation) and death events. The event that actually occurs corresponds to the minimum of 
these durations, over the set of all individuals. 



3.2 Examples 

Simulations of particular cases are now presented. For each example, simulations have been run several 
times ; only one particular simulation is presented for each example as an illustration. The programs 
have been written with the R freeware (0) and Matlab. 



3.2.1 Example 1: Evolution of offspring size in age structured populations with size- 
dependent competition 

In this example, we are interested in a population with size-dependent competition, and ask how the 
evolution of body size changes when we take the growth of individuals into account. 

Individuals are characterized by their body size at birth xq € [0, 4] which is the heritable trait subject 
to mutation, and by their physical age a € [0, 2], with aging velocity 1. Body size is an increasing function 
of age: 

X — xq + g a, (20) 

where g is the growth rate, which is assumed constant and identical for all individuals. 
An individual of trait xq e [0, 4] gives birth at the age independent rate: 



b{xo) 



Xq. 



(21) 



This is a simple expression of the standard trade-off between fecundity and offspring size {e.g. Stearns 
[82|). With probability p = 0.03 a mutation occurs and affects xq. The new trait is Xq — min(max(0, xq -\- 
2/2,q),4), where is a Gaussian r.v. with expectation and variance 0.01. With probability 1 — p, the 
offspring inherits its parent's trait, xq. 

The death rate of an individual with trait xq G [0, 4] and age a e [0, 2] living in a population 
Z e A^p([0,4] X [0,2]) is given by: 



d{xo, a,Z) = IU{xo + ga 
Jx 



x'q — g a)Z{dx'Q, da), 



(22) 



where 



U{x-y) =- 



300 



1 - 



1 



1 + 1.2 cxp(-4(x- ?/))^ 





2 




300 



(23) 



is the asymmetric competition kernel introduced by Kisdi [55, which gives a competitive advantages to 
larger individuals. The important point is that the size of competitors, hence the intensity of competition 
between them, are not constant when g 0. Even in a monomorphic population, size varies across 
individuals and competition experienced by each individual varies with its size, which is age-dependent. 

In the simulations of Figure [21 we chose g = 0, g = 0.3 and g = 1. The population at time t — 
is monomorphic with trait xq — 1,06, and contains Nq — 900 individuals. The initial age distribution 



^http:/ /www. r-project.org/ 
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Figure 2: Trait dynamics in an age structured population with size- dependent competition. The evolving trait is 
offspring size. See Example 1 in main text for details, (a) g — 0, (b) g — 0.3, (c) g = 1. 

is uniform on [0, 2], and the corresponding sizes are computed using (|20p . For each t the support of the 

measure J2i=i ^(xo)i represented. 

For g = (constant size during hfe), we observe a branching phenomenon around t = 300: the distri- 
bution of trait xq sphts into two "branches". The initial branch first stabiHzes around some equihbrium 
[xq ~ 2.7), and then, a second branch appears below [xq ~ 2). For g = 0.3, the population branches into 
two subpopulations {xq ~ 2.7 and xq ~ 2) sooner, around t = 100. For <? = 1, the branching phenomenon 
is not observed anymore. The occurrence of branching in Figl^] (a) shows that the coexistence of subpop- 
ulations can be promoted by the competition-fecundity trade-off: small individuals may compensate for 
their competitive inferiority by reproducing more often. In (b), non-zero body growth results in a wider 
size-structure in the population, so that the conditions of branching are met earlier. Here, growth acts 
as a mixing factor, which lessens the differences between traits. In (c), individual growth is so fast that 
the differentiation of a subpopulation of smaller individuals is always prevented by strong competition 
with very large individuals. 

3.2.2 Example 2: Coevolution of offspring size, growth rate and age at maturity under 
size and stage-dependent competition 

As in Example 1, we consider an age- and trait-structured population where the traits are the size at 
birth xq £ [0, 4], the growth rate g € [0, 2], and the age at maturity gm G [0, 2]. The scalar age a e [0, 2] 
(with aging velocity 1) still corresponds to the physical age. 

The age at maturity om separates the life history of individuals into two periods: a growth period 
and a reproduction period. Calsina and Cuadrado [l3j also studied a population that is structured into 
two age classes, juveniles and adults. Assuming that the length of the juvenile period is exponentially 
distributed, they described population dynamics by a system of two ODEs. They then modelled the 
evolution of the age at maturity by considering whether a small mutant population can invade the 
resident population, leading to questions of stability of the stationary solutions. Ernande et al. [34| 
considered a similar problem: their population was divided into age classes, and the age at maturity was 
plastic and determined by a heritable reaction norm that was subject to mutation. 

In our model, the size of an individual with trait {xq, g, qm) and age a is given by: 

x{a,xo,g,aM) ^ xo+ g[a haM)- (24) 

Before maturity, individuals (called juveniles) do not reproduce and invest all their resources into growth. 
Once the age ajv/ has been reached, individuals (called adults) have constant size and begin to reproduce 
(see FigEl). 

The birth rate of an individual with traits (xcg, ajv/) G [0,4] x [0,2] x [0,2] is for juveniles and a 
decreasing function of xq for adults: 

6(a;o,aM,a) = (4- xo) la>aj,f (25) 
The traits xq, g and aM are heritable and can be affected in offspring by mutation that occur with 
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Figure 3: Size x as a function of age a. The age at maturity is denoted by um- Before maturity, juveniles born 
with size xo grow at rate g. After maturity, the size of adults is constant. 

probability p — 0.03. Mutated traits are chosen according to: 

= min(max(0, + yxa)A) s' = niin(max(0, g + Vg),2), a'^ = min(max(0, hm + 2), 
where Uxq, Vg and are Gaussian independent r.v. with expectation equal to and variance equal to 



The death rate of an individual with traits {xo,g, anj) G [0, 4] x [0, 2] x [0, 2] and age a e [0, 2] in the 
population Z e A4p([0, 4] x [0, 2] x [0, 2] x [0, 2]) is given by: 



{do X 5 + 10 ^ ^^U{xq ^- g a — x'q — g' [a! /\ a'j^j))Z{dx'Q,dg' .da'j^i^da') if a < om, 
[/(xo + goA/ - - g'{a' A a'j^j))Z{dx'Q,dg' ,da\,f,da') if aM < a < 2, 
+00 if a > 2, 



where U is the competition kernel introduced in (j23p and da > measures the adversity of environmental 
conditions. During the juvenile period the main factor of death is related to growth: individuals with high 
growth rates g need more ressources and incur a survival cost of foraging ; thus, the density- independent 
component dog of their death rate is higher if do is larger. Then, during the adult period, there is no more 
growth and death is fully determined by logistic competition, which was soften for juvenile individual. 
This may apply to species that shift to different habitats or resources as individuals reach maturity. 

As in Example 1, we chose as initial state a monomorphic population {xq = 1.06, g = 0.74 and 
ajv/ = 0.20) with A^o — 900 individuals whose age is uniformly drawn in [0, 2]. Fig. [4]shows that changing 
do can result in very different dynamics. When do is small (Figl^a), do = 3.310^"^), individuals with 
higher growth rates do not pay much of a cost (Vg € [0, 2] , < dog < 6.6 10~^) . The juvenile death rate is 
small (of the order of 10~^ for a population of 1000 individuals) compared to the adult death rate (of the 
order of 1 ) . Since competition favors larger individuals (the term U{xo+gaM — x) describing competition 
exerted by an individual with size x on an adult with size xq + gaM is smaller when xq + gaM is large) . 
Hence there is a strong selective advantage to grow fast {g ~ 2) and long (om — 1-3). The reproduction 
loss caused by delayed maturity is further compensated by the evolution of relatively small size at birth. 
Consequences for population dynamics are shown in Fig[Sl The population stabilizes numerically around 
a state with high growth rate and late maturation. 

When do is large (FigHl^b), do = 0.5), a state with moderate growth and early maturation evolves. 
We observe that the growth rate decreases towards 1 and evolution seeks to minimize the age at maturity 
(om — 0.3). To understand the changes when do increases, let us introduce the probability that an 
individual of traits (xo,(7, om) born at t = c survives until age at maturity um in a population (Zt)tgR^: 



When the population size is bounded by N, the probability n(a;o, g, um) is upper and lower bounded by: 



0.01. 




(26) 




(27) 



(a) 



(b) 




Figure 4: Time dynamics of three co-evolving traits: size at birth, growth rate, age at maturity (rows 1-3), and 
population structure (distribution of maternal ages at birth) in the "equilibrium state" (row 4)- See Exemple 2 in 
main text for details, (a): do = 3.3 10~^, (b): do = 0.5. 



In the simulations, we observe that the size of the population remains bounded and that it stabilizes 
around an "equilibrium" value (600 for do = 3.3 10^^ and 200 for do = 0.5, see Fig. [5]). The "equilibrium" 
traits obtained in the two simulations of Figure [3] are distributed around [xq, g, um) = (1-8,2, 1.3) and 
{xo,g,aM) = (2.6,0.8,0.3). The strategy that is observed can be naturally explained. For do = 0.5 
and N = 600, we hence obtain 0.27251 < n(1.8, 2, 1.3) < 0.27254, instead of n(1.8, 2, 1.3) ~ 0.99 when 
do = 3,3.10-3, whereas 0.88690 < 0(2.6,0.8,0.3) < 0.88692. For do = 0.5, the traits {xo,g,aM) = 
(2.6, 0.8, 0.3) are more competitive. They correspond to individual who have a higher probability to 
reach maturity and reproduce, even if these individuals arc likely to be exposed to higher competition 
pressure that individuals with trait (2:0,5, cm) = (1.8,2,1.3). 

In the case of a favorable environment with abundant ressources (do is small), the growth period has a 
small cost and a long growth phase (om large) allows individuals to reach a mortality refuge: i.e. a state 
where they escape the effects of strong competition pressure. When dp increases, the mortality refuge 
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becomes inaccessible: long period of initial growth result in low probabilities of survival until reproduction. 
An alternative life profile appears, with early reproduction (small um) and small individuals. Similar 
phenomena have been discovered and discussed in Taborsky et al. [83^ . 




(a) (b) 

Figure 6: Long term density of maternal age at birth (a) and lifespan (b) with respect to the mortality cost of 
growth, do. Horizontal projection contours show that as do increases, the distribution of maternal age at birth 
changes smoothly (a), whereas the distribution of lifespan undergoes an abrupt transition (around do = 0.2) (b). 
See main text for details. 



The transition from "cheap" growth (low do, Fig- SJa)) to "expensive" growth (high doj Fig. [IJb)) 
is investigated numerically in Fig. [51 where c?q varies between 0.01 and 0.6. Figure [S]^a) shows the 
kernel density estimator [e.g. 8J) of the maternal age at birth. Fig. [Htb) shows the density estimator 
of lifespan. When do increases from 0.01 (favorable environment) to 0.6 (hostile environment), lifespan 
and reproduction age decrease continuously from 1.5 to 0.5. Compared to the results of [HS], the surfaces 
displayed in Fig. [6] do not show strict discontinuities, but the smooth variation of the distribution of 
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maternal age at birth (a) contrasts with an abundant shift in the distribution of hfespan (b). This shift 
occurs in spite of the survival probability ()26p (and hence the lifespan density) being a continuous function 
of do. 

4 Conclusion 

Classical models for the dynamics of populations with continuous age structure ignore stochastic processes 
operating at the level of individuals, and therefore take the form of deterministic PDEs, e.g. the celebrated 
McKendrick-von Foerster equation. Stochastic models that start at the level of individual processes and 
account for interactions between individuals face serious mathematical challenges. One way forward is to 
derive large population limits from these IBMs. This approach was followed here for a large class of IBMs, 
for which the large population approximation is a PDE that generalizes the Mc Kendrick-Von Foerster 
model. There are qualitative limitations of the large population approximation, that we demonstrated by 
showing almost sure extinction of the logistic age-structured population process. However, stochastic and 
deterministic approaches appear to yield complementary insights into the population dynamics: prior to 
extinction, stochastic trajectories are proved to spend an exponentially distributed time near the solution 
of the deterministic approximation model. In general, the rigorous derivation of the deterministic limit 
from stochastic individual-level processes allows one to derive confidence intervals for model parameters 
and thus open new ways of confronting models to real individual data. 

In populations that are structured by age and trait variation (due e.g. to genetic mutation), the 
analytical study even of the large population approximation becomes intractable. But the rigorous con- 
struction of the individual-based stochastic process yields an efficient algorithm for numerical simulations 
of population age and trait distributions. Two biological examples were presented. In the first example, 
offspring size varies genetically and evolves under size-dependent competition and the genetic constraint 
of a tradeoff with the birth rate. Individual growth and the resulting age and size structures have dra- 
matic influences on trait evolution. Moderate growth shapes the population size structure in a way that 
exacerbates competition and favors the rapid split of the population into two 'evolutionary branches' 
{i.e. the trait distribution becomes bimodal). With more rapid individual growth, the size distribution 
widens even more and competition intensifies to the point where the divergence of trait branches becomes 
impossible {i.e. the population trait distribution remains unimodal). 

In the second example, three traits were allowed to vary genetically and co-evolve: offspring size, 
growth rate, and age at maturity. The corresponding population model assumes size- and stage-dependent 
competition. As in the previous example, offspring size evolves under the genetic constraint of a trade-off 
with the birth rate. The growth rate evolves under a risk-competition tradeoff : faster growth (requiring 
e.g. acquisition of more resources, hence riskier behavior) entails a mortality cost but ensures a larger 
size at maturity and hence a competitive advantage in the reproductive stage. Age at maturity evolves 
under a reproduction-competition tradeoff: reproductive {i.e. mature) individuals face a higher mortality 
risk due to competition than juveniles. Model simulations were performed to analyze the effect of the 
mortality cost of growth on the traits' coevolution. High cost promotes the evolution of a moderate 
growth rate, together with large size at birth and very early maturity. The corresponding distribution of 
longevity in the population is skewed towards low values. Interestingly, as the mortality cost of growth 
decreases, the longevity distribution shows a relatively abrupt shift towards large values - a qualitative 
change in the population demography that was not anticipated given the continuous dependence of the 
longevity density on the cost of growth. 

These examples highlight that the mathematical and numerical analysis of stochastic population 
models with age and trait structure have the potential to uncover unexpected phenomena of biological 
interest and to advance ecological and evolutionary population theory in significant ways. 
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